Stark effect and generalized Bloch-Siegert shift in a strongly driven two-level system 
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A superconducting qubit was driven in an ultrastrong fashion by an oscillatory microwave field, 
which was created by coupling via the nonlinear Josephson energy. The observed Stark shifts of 
the 'atomic' levels are so pronounced that corrections even beyond the lowest-order Bloch-Siegert 
shift are needed to properly explain the measurements. The quasienergies of the dressed two-level 
system were probed by resonant absorption via a cavity, and the results are in agreement with a 
calculation based on the Floquet approach. 
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The characteristics of matter and light become inter- 
twined upon interaction. The interconnection of the two 
can be observed in atomic and optical physics by set- 
ting atoms inside mirrored cavity resonators, whereby 
coupling the zero-point vibrations of the field to atomic 
transitions. Pioneering studies of quantum physics have 
been carried out in the scheme (see, e. g., [1]). 

Recently, an increasing emphasis has been put on the 
study of analogous physics in a setting of electromagnetic 
modes interacting with discrete systems, but in a solid- 
state environment. Such quantum few-level systems, or, 
artificial atoms, have been implemented based on quan- 
tum dots and superconducting Josephson qubits [2-4]. 

Apart from zero-point vibrations in a cavity, natural 
or artificial atoms can be coupled to a driven laser field. 
One of the effects of the field is the dynamic (ac) Stark 
shift [5-9] of the energy levels. For off-resonant driving, 
the shift scales linearly in the number of quanta in the 
field. An additional correction, the Bloch-Siegert shift, 
appears for an oscillating, rather than circularly polar- 
ized field [10, 11]. A strong drive with the Rabi frequency 
Un becoming a sizable fraction of the atomic level spac- 
ing is needed to reveal the realm of such delicate phe- 
nomena, which poses a challenge for experimenting with 
real atoms. A good understanding of the physics is im- 
portant also in the sense that the artificial systems have 
been actively investigated due to their promise of setting 
up quantum information processing. 

In the present work, we investigate the energy levels 
of an artificial two-level system driven by an oscillatory 
field originating from a harmonic drive via the Joseph- 
son energy. For the purpose, we have developed a qubit- 
resonator setup, where matrix elements allow to carry out 
the measurement over the entire excursion covered by the 
drive. We measured large Stark shifts of the qubit level 
spacing unseen in atomic systems. The shifts are found to 
exhibit unconventional and to some extent nonmonotonic 
dependence on the field amplitude. This work is the first 
observation of Bloch-Siegert type of correction in driven 
systems other than atomic systems [12] or magnetic reso- 
nance [11]. Unlike all but few previous experiments [13], 
we have to go beyond the usual lowest-order such cor- 
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FIG. 1: Schematic of the experimental setup, (a) The single- 
Cooper-pair transistor qubit consists of a superconducting 
loop interrupted by two Josephson junctions that separate 
a small island. The island has two relevant charge states cor- 
responding to n or n + 1 Cooper pairs; (b) Except near the an- 
ticrossings, the energy eigenvalues depend nearly sinusoidally 
on the applied flux $. The flux is driven likewise sinusoidally 
in time; (c) A more detailed view of the circuit shows the 
microwave reflectometry readout via an LC-resonator formed 
by on-chip lumped circuit elements [27]. 



rection in order to account for the data. This work dif- 
fers from the previous studies on strong drive with su- 
perconducting or quantum dot qubits, which focused on 
the Landau-Zcner-Stiickelbcrg (LZS) effect [14, 15], and 
where the coupling of the drive to the qubit Hamiltonian 
was linear [5, G, 16-26]. 

We use as the two-level system a single- Cooper-pair 
transistor (SCPT) [2, '-i], a superconducting qubit which 
consists of two small-area tunnel junctions with the 
Josephson energies Eji and Ej2, see Fig. 1. The phase 
difference across the SCPT can be tuned by means of a 
magnetic flux $ applied through a superconducting loop. 
In order to minimize background charge noise, the gate 
voltage is adjusted so that two charge states differing by 
one Cooper pair are degenerate [3]. Neglecting higher 
charge states, the Hamiltonian of the SCPT written in 
the charge basis consists of the Josephson coupling ener- 
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The parameters are the total Josephson energy Ejo ~ 
Eji + Ej2, the asymmetry d = {Eji - Ej2)/Ejo, and 
the flux quantum $o = h/2e. The energy eigenvalues 



Ej{^) = £;jo^cos2(7r$/$o) + d2sin^(7r$/$o) of ff, (1) 
are depicted in Fig. lb. 

The applied flux $ = $b + $l cos{ujLt) consists of a 
static bias $h and of a time-dependent part with the am- 
plitude $L which is analogous to the intense laser field 
in atomic physics. The on-chip fiux coil was designed to 
have a large mutual inductance of ^ 5 pH to the super- 
conducting loop in order to achieve a desired fast control 
of the flux bias over the span of several $o- Owing to the 
nested sinusoidal time dependence of the energy, as well 
as unequal coupling to the different cr^^y^z, the driving 
field substantially deviates from a circular polarization, 
and novel phenomena appear. 

In the following we find the coupled states of the qubit 
and the field. This leads to the picture of the dressed 
states formed by the qubit and the light field. The spec- 
trum is obtained as quasienergies, which can be consid- 
ered as the characteristic energies of the combined system 
of the qubit and the field [1]. They repeat periodically 
at intervals lul as depicted in Fig. 2. 

Assuming d small, we diagonalize the first term in 
Hq (1), replacing a^- — >■ o-^, ay -^ ax- The eigenstates 
of this term are called "bare" below. For simplic- 
ity, we give explicit formulas in the case where we 
keep only the 0th and 1st harmonics of the driven 
Hamiltonian (1). The numerical calculations, however, 
are done for the full Hamiltonian. To account for 
the temporal variation of the longitudinal field, we 
transform to a rotating frame. The transformation 
is obtained hy U — exp ( — iazAsm{ujLt)/2huj]^), 
where A = 2£',7o Ji(7r$L/$o) sin(7r$b/$o)- We 

define hujQ = -£'joJo(7r$L/$o) cos(7r$b/$o), 

B = £;jodJo(7r$L/*o)sin(7r$b/$o) and fife = 
[B + kujLdcot [Tr^b/^o)]Jk{A/'hujL) and find 
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The strongest effect of the "laser" field is a mod- 
ification of the qubit splitting to Tiuiq by a factor of 
Jo(7r$L/$o)- This can be understood as a rectification of 
the drive by the qubit nonlinearity which effectively shifts 
the bias point, as seen by the transition from the dotted 
to dashed lines in Fig. 2. The positions of the resonances 
(level crossings), hence become shifted, cf. Fig. 2. 

Apart from that, the fife terms describe transverse cou- 
plings induced by the time-dependent laser field. To take 
that into account, we consider _ffid (2) as a matrix oper- 
ating on the longitudinally dressed states jcr, n) [IG, 21]. 
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FIG. 2: Calculated quasienergies as a function of the bias 
flux $5 with the experimental parameters. The fully dressed 
states (red lines) are the eigenvalues of the Floquet matrix. 
The quasienergies repeat periodically in energy with period 
TiujL- The level splittings are denoted by ?iA and hiujL ~ A). 
The vertical green and blue arrows denote transitions induced 
by 3.5 GHz probe field. We also show the quasienergies of the 
bare states (c&l = 0, dotted lines) and longitudinally dressed 
states (rifc = 0, dashed lines). 



They are labeled by a qubit ('spin') index a =tii and 
a Floquet ('photon') index n. The couplings have an 
effect of opening a gap at the fcth anticrossing. This 
shifts the quasienergies from their uncoupled values (the 
ac Stark shift [7]). For very weak couplings, Hfe <^ ujl, 
it is enough to keep only the couplings between resonant 
levels in the Floquet picture (rotating-wave approxima- 
tion, RWA). RWA allows for analytic solutions and gives 
the Stark shift of ywSl+TIJ — \6l\, where 6l is the de- 
tuning from resonance. Among other approaches, RWA 
was sufficient for description of recent LZS experiments 
on superconducting qubits [10, 18, 21]. 

In our case, Q.k/^L is closer to 1, depending on the 
dc and ac bias, and thus a few orders of the perturba- 
tive expansion in this parameter beyond RWA should be 
accounted for (the second order, similar to the Bloch- 
Siegert contribution, and beyond). Adding to i?id (2) 
the photon energies, we solve numerically for the eigen- 
values and eigenstates of the relevant part of the infinite 
Floquet matrix [2n]. The calculated quasienergy split- 
ting A for the drive amplitude $l = 0.26 $0 is plotted 
in Fig. 2. Over the full ($b,<i>L) parameter space, the 
splitting is shown by contours in Fig. 3(b). 

In Fig. 3(b) the Rabi resonances appear as light and 
dark tracks, corresponding to A/wi w 1 or 0, respec- 
tively. These resonances start as vertical lines at small 
<i>i, but curve to the right with increasing $i. The 
curving is a consequence of approaching the first zero 
of Jo(7r<I>L/<I>o) at $i w 0.77$o, which implies a dy- 
namic collapse of the qubit. This effect of the non- 
linear longitudinal driving bears similarity to the co- 
herent destruction of tunneling [29] by transverse driv- 
ing. It is visible as the dark horizontal band around 
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FIG. 3: (a) The measured resonance absorption of the probe 
signal plotted in the $(, — $i, plane. Dark corresponds to 
higher absorption. The vertical scale is obtained by compar- 
ison to theory; (b) The landscape of the quasienergy split- 
ting A calculated from the Floquet matrix of Hamiltonian 
(2) with the experimental parameters. The fc-photon Rabi 
resonances appear as light (odd k) or dark (even k) tracks 
that are marked by index k at small $i and curve to the 
right with increasing i1>_l. Fig. 2 shows a cut along the ho- 
risontal line denoted by the black arrow. The blue and green 
lines denote the resonances in Eq. (4). The solid dots are the 
experimental resonances picked up from panel (a), or from 
the measured phase shift (two rightmost lines). 



^L ~ 0.75$o in Fig. 3(b). The plot also shows sev- 
eral conical points (crosses) , where the energy values are 
degenerate (A/wl = 1 or exactly) [30]. 

The experiments were performed in a dilution cryostat 
at a temperature of 30 mK. We first measured the sam- 
ple parameters independently on the strong-drive experi- 
ments. Especially near the charge-flux degeneracy point, 
detailed mapping of the energy landscape is required in 
order to later obtain accurate comparison to the expected 
frequency shifts. The values of the total Josephson en- 
ergy Ejo/h = 27.0 GHz and d = 0.19 were obtained by 
microwave spectroscopy at a very small drive, which di- 
rectly yields the bare qubit's level spacing by making a flt 
to numerically evaluated energies taking into account 10 
charge states. The ratio Ejq/Ec = 8.0 was obtained as a 
fit to the ground-state inductance over the V^-$b plane. 
Here, Ec ~ e^/2(Ci + C2), and the parallel junction ca- 
pacitances are Ci and C2. The ground-state response 
agreed with the spectroscopic measurement of Ejq and 



Instead of measuring directly the excited state popu- 
lation [18], we probe the quasienergies with a weak fiw- 
signal. The probe is produced by coupling the qubit with 
an LC-resonator via the total flux $ in the loop (Fig. Ic), 
and the information is encoded in the exchange of energy 
with the qubit, or in the dispersive frequency shift. The 
resonator is formed by the inductance of the supercon- 
ducting loop, L ~ 410 pH, and by the lumped element 
capacitors C ^ 10 pF made out of Al oxide between the 
Al films. Since the LC frequency ~ 3.5 GHz is smaller 
than the minimum level spacing dEjo r^ 5.1 GHz of the 
bare qubit, the qubit-resonator interaction does not no- 
ticeably influence the energies. 

For measurement, the circuit is excited at a frequency 
close to the resonance frequency, and the phase and am- 
plitude of the reflected signal are recorded. We devel- 
oped the following semiclassical description of the mea- 
surement. The resonator coupling is equivalent to adding 
a resonator-induced flux $p(i) = <I>p cos(a;pt) into $ = 
<&6 + <&L(i) + $p(0- We assume that ^p <C $0 so that it 
does not perturb the dressed states significantly. Trans- 
forming to the basis of longitudinally dressed states, the 
Hamiltonian is H = Hi^ + Hp, where Hi^ is given in (2), 
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We have defined 5n — Tid (wq + nujL)Jn{A/hujL). Assum- 
ing small $p, the probe Hamiltonian (3) can be treated 
as a small perturbation. Besides the pcrturbative analy- 
sis, we have run a full simulation on the rcfiection mea- 
surement. For this purpose, we have derived the equa- 
tions of motion for our circuit using quantum network 
theory [31, 32]. The resonator is treated classically, and 
it is coupled to Bloch equations describing the qubit. 

In the dispersive regime, where the resonator is sub- 
stantially detuned from the qubit, the resonator fre- 
quency changes according to the Josephson inductance 
of the qubit, which was used to obtain the parameters of 
the qubit in the ground state. With strong drive, how- 
ever, the dispersive signal becomes overwhelmed by the 
absorptive response. Because the probe frequency up is 
smaller than the laser frequency lol (wL/27r = 6.11 GHz 
and LOpjI'K = 3.5 GHz), we can see two different transi- 
tions in the absorptive measurement. These are the clos- 
est states, corresponding to energy differences ?iA and 
'hiujL — A) in Fig. 2a. The resonance conditions are 



Up (blue) 



WL - Up (green) , (4) 



where the color coding is related to Figs. 2 and 3a. They 
can be interpreted as the lowest transition in a fluorescent 
triplet, and the transition between Rabi states [1]. 

In the $b-$L plane, the resonance conditions appear 
as contours (Fig. 3b) . The matrix elements for the tran- 
sitions are determined by Up (3). Since they are nonzero 
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level spacing. In order to obtain a good agreement be- 
tween experiment and theory, we have to consider higher 
orders beyond the rotating-wave approximation which 
amounts to including a generalized Bloch-Siegert shift. 
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FIG. 4: The shifts of the spectral lines due to the qubit 
driven by a strong field. The data points are from Fig. 3. 
The theoretical curves were produced with the rotating wave 
approximation in the adiabatic basis (red lines), or with 
the full numerical result (green). In (c), the dashed line is 
the analytical solution for the resonance condition, {kujL + 
ajp){[Jo(7r$i/$o)]"'-l}, fc = 3. 



almost everywhere, accurate mapping of the energy land- 
scape is possible. The resonant energy flow [27] from the 
resonator to the qubit can be seen as increased absorp- 
tion as well as a phase shift in the reflection measurement. 
The former is shown in Fig. 3a. The measured locations 
of the absorption maxima follow closely the resonance 
lines (4) as seen by overlaying them with the theory in 
Fig. 3b. The full simulation of the measurement is in 
good agreement with the measured reflection. 

One can also extract from Fig. 3 the actual Stark shift 
of the spectral line. The shift from the undriven case is 
illustrated in Fig. 4 for the three resonances labeled A-C 
in Fig. 3. The overall trend is that the shift grows with 
the drive amplitude as [Jo(7r$L/<i>o)]~ — 1 due to the 
rectification, which is obtained by inverting the effective 
qubit splitting wq for comparing to the original splitting. 
Moreover, the 'Rabi' gap ilk modulates the shift near 
the kth resonance. Its Bessel-type dependence on the 
drive gives rise to the nonmonotonic Stark shift in Fig. 4. 
For the curve C the relevant coupling fl^ is weak near 
^b/^o = 1 and can be neglected. It is hence enough to 
account for the rectified level splitting ujq in this region, 
cf. Fig. 4c. 

The data display a good agreement to the full theory 
calculations (green). We also show a comparison to the 
RWA calculation (red), which clearly fails to explain the 
data. For the resonance of Fig. 4a, higher order correc- 
tions even up to fourth order still substantially deviate 
from the data. The difference between the full calcula- 
tion and the RWA can be considered a generalized Bloch- 
Siegert shift. 

To conclude, we have measured and theoretically an- 
alyzed a superconducting qubit strongly driven via the 
Josephson energy. The spectral lines of the qubit expe- 
rience a remarkably strong Stark shift up to the original 
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